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Abstract 

We consider a perfectly homogeneous and isotropic universe which undergoes a sudden 
phase transition. If the transition produces topological defects, which we assume, pertur- 
bations in the geometry and the cosmic fluid also suddenly appear. We apply the standard 
general relativistic junction conditions to match the pre- and post- transition eras and 
thus set the initial conditions for the perturbations. We solve their evolution equations 
analytically in the case when the defects act as a coherent source and their density scales 
like the background density. We show that isocurvature as well as adiabatic perturbations 
are created, in a ratio which is independent of the detailed properties of the defects. We 
compare our result to the initial conditions currently used in the literature and show how 
the cosmic fluid naturally "compensates" for the presence of the defects. 
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1. Introduction 



The question of the generation of the small cosmological inhomogeneities at the origin 
of the large scale structures observed in the universe is still open. Two main approaches, 
each with its troop of specific models, are currently in competition. On one hand the infla- 
tionary scenario (see e.g. [1]) explains cosmological inhomogeneities by the amplification, 
due to accelerated cosmic expansion, of inescapable quantum fluctuations. On the other 
hand, the topological defect scenario, which is based on the idea of spatially differentiated 
spontaneous symmetry breaking (see e.g. [2]), explains these inhomogeneities by the ap- 
pearance of topological defects which drive fluctuations in other types of matter. Future 
observations of the small scale anisotropies of the CMBR (Cosmic Microwave Background 
Radiation), in particular by the planned MAP and PLANCK satellite missions, should 
discriminate between the two scenarios. 

In the present work, we shall be interested in the second one. Whereas the calculation 
of cosmological perturbations in inflationary models, at least the simplest ones, seems 
to have now reached a stage of maturity and clarity, there is still some confusion in the 
literature on perturbations seeded by topological defects, even on the question of how to 
set the initial conditions. By initial conditions, we mean, as is usual in cosmology, the 
state of cosmological perturbations for the various matter species at a past epoch during 
the radiation era when all scales of cosmological interest today were larger than the Hubble 
radius. Although far in our past, this initial epoch is taken in general very long after the 
phase transition supposed to have given birth to the topological defects. The reason is that 
it is difficult to trace numerically the evolution of perturbations on a very long duration. 
The problem is thus to translate the "starting" conditions imposed by the phase transition 
into "initial" conditions at the time the numerical computations begin. 

At this "initial" time, the topological defects are thus supposed to already exist. The 
main difficulty is then to set the initial value of the perturbations of the cosmological fluids 
which arc compatible with the distribution of the defects. Several methods have been 
dcvclopped to determine and then implement these compatibility conditions: "integral 
constraints", "compensation", "pseudotensor" . . . Unfortunately, these methods are rather 
intricate and confusing... 

The purpose of the present paper is to present a self-contained, and hopefully clear, 
analytical derivation for the setting of these initial conditions. To do so, we start our 
analysis before the phase transition, when the universe is remarkably simple since it is 
supposed to be strictly homogeneous and isotropic. Starting from this pre-transition state, 
one should, in principle, study the detailed evolution of the scalar field at the origin of the 
defects. An important simplification arises from the fact that the defects can be considered 
as perturbations and their evolution be obtained by solving their equations of motion in 
the background geometry. However the problem remains very complicated and requires 
heavy numerics. What will be retained for our purpose is that the quantities describing the 
defects can be seen as "external" sources for the evolution equations of the perturbations 
of the cosmic fluids. 

In order then to relate the pre-transition Friedmann Lemaitre Robertson Walker 
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(FLRW) spacetime to the subsequent perturbed universe containing defects, we shall as- 
sume a sudden phase transition where the defects are instantaneously "turned on". Ein- 
stein equations then impose some matching conditions between the perfectly homogeneous 
spacetime and the defect populated one. These matching conditions provide constraints 
but are not sufficient in themselves to determine completely the subsequent state of the 
perturbations. 

If one wishes to say more about the post-transition perturbations, one must make some 
assumptions about the defects. The emergence of defects can be seen as a random process 
and the information on defects will thus be of a statistical nature. We shall assume that the 
sources satisfy the properties of causality, scaling and coherence. Causality is physically 
required and simply states that the sources are uncorrelated on scales larger than the 
Hubble radius. Scaling means that the statistical properties are invariant in time up to a 
rescaling with respect to the Hubble radius. This property is a simplifying one but can be 
justified by the convergence towards scaling observed in numerical simulations. Here, we 
assume scaling immediately (that is within one Hubble time) after the phase transition. 
Finally, coherence is a very stringent assumption, because it means that all the statistical 
properties of the sources are reduced to their statistics at a given time, the time evolution 
being completely deterministic. However, to consider coherent sources is more general than 
it seems because any source for our purpose can be decomposed into a sum of coherent 
ones. 

The advantage of considering coherent sources is that, knowing the time dependence 
of their correlators from scaling and causality, one can then solve for the evolution of the 
perturbations by a simple use of Einstein's equations. In the long wavelength limit, i.e. 
for scales larger than the Hubble radius, one obtains the solution explicitely as a sum of 
power-law terms. 

From our explicit analytical solutions, we are able to justify some of the statements 
which can be found in the literature. In particular, the solutions for the perturbations, 
which can be decomposed into terms driven by the sources and "homogeneous" (in the sense 
of differential equations) terms, are shown to be dominated by the source driven terms. 
We also recover on our solutions the so-called phenomenon of "compensation". Another 
conclusion of our work, which is original to our knowledge, is that, when one allows for cold 
dark matter, the sources will drive the long-wavelength perturbations into a combination 
of adiabatic and isocurvature perturbations, the relative ratio being a universal constant. 

The paper is organized as follows. In section 2, we write the linearised Einstein 
equations. Our formalism is based on the usual cosmological perturbation theory using 
Bardeen type gauge-invariant quantities, where one distinguishes between scalar, vector 
and tensor perturbations. In addition to the usual geometrical and perfect fluid type 
matter variables, one must introduce the energy momentum tensor of the defects. This 
tensor is automatically gauge invariant because, as stated before, the defects are considered 
as perturbations. In Section 3, we describe the phase transition. In Section 4 we review 
the matching conditions on a constant energy density surface and apply them to the case 
of a sudden phase transition in the universe. They enable us to make the link between the 
unperturbed universe and the post-phase transition perturbed universe, or rather they give 
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constraints that physically admissible configurations must satisfy. Section 5 deals with the 
statistical properties of the sources and the notions of coherence, scaling and causality are 
introduced in detail. In section 6, we give the behaviour of perturbations larger than the 
Hubble radius after the transition, which will constitute the "initial" conditions. Finally, 
section 7 comments on the obtained solutions. 

2. Linearised Einstein equations 

2.1 The background 

The universe at large appears to be remarkably homogeneous and isotropic and gov- 
erned by the gravitational force created by its material content, to wit a mixture of ra- 
diation and dust. It is therefore well described by a Robertson- Walker geometry whose 
time evolution satisfies Friedmann's equations. Phase transitions [3] [2] can have occured 
in the very early universe when it was pure radiation and its spatial curvature negligible. 
We thus take the background line element in that era to be : 

ds^ = a^{r]){-dr]^ + Sijdx'dx^) (1) 

where = t] is conformal time, i = 1,2,3, three cartesian coordinates, and a{r]) the 
scale factor. The Friedmann equations are : 

n' = -n^ , npa^ = STt^ (2) 

where a prime denotes a derivative with respect to conformal time, where Ti. = a' /a, where 
p is the radiation energy density and n = SttG is Einstein's constant. The solution of (2) 
is : 

.,1 2 3 p' 4 

n = - ^ a = aorj , Kpa = — = — (3) 

7] rj'' p rj 

where oq is an integration constant which can be chosen so that a = 1 today. 

2.2 Scalar perturbations 

Following Bardeen [4] we split the perturbations of the geometry and the matter 
variables into "scalar" , "vector" and "tensor" parts. In this paragraph we define the scalar 
part and write the linearised Einstein equations for the gauge invariant scalar perturbations 
(for reviews of this formalism, see e.g. [5-7]). 

The line element of a perturbed Robertson- Walker space time reads, when the per- 
turbations are scalar and the background given by (1-3) : 

ds^ = alrf [-(1 + 2A)drf + 2diBdx'dr} + {(1 + 2C)5ij + 2dijE} dx'dx^] (4) 

where A, B, C, E are four "small" functions of space and time. Under an infinitesimal 
coordinate transformation 

r]^r] + T , x''^x'' + d''L (5) 
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where T and L are two arbitrary first order functions of rj and a;*, the four scalar metric 
perturbations A, B, C, E transform as 

A + T' + -, B^B-T + L', C^C+-, E^E + L. (6) 
V V 

One can thus introduce two gauge independent scalar metric perturbations, for example : 

^ ^ -C --(B - E') , <^^A+-(B-E') + (B-Ey. (7) 
V V 



The energy-momentum tensor of the matter content of this perturbed universe can be 
written as : 

+ ST^, + e^,. (8) 

T^i, {jJL^v = 0,1,2,3) is the energy-momentum tensor of the homogeneous and isotropic 
radiation background ; ^T^^y is its perturbation : its scalar components can be expressed 
in terms of two scalar quantities, S = Sp/p, the density contrast, and v, the velocity 
perturbation, as (see e.g. [5-7]) : 

^{6 + 2A) , n6Ti = -^dJB+i 
T]^ \ 3 



kST^o = —{S + 2A) , nSTi = -^di(B+^v 

"l " ' 

nSTf;^ = — [6,j{2C + S) + 2d,,E] . 



(Note that we describe the material content of the universe as a single radiation fluid (with 
no anisotropic stresses). This is justified since at the era of the phase transition all matter 
is highly relativistic.) 

In the coordinate transformation (5) 5 and v transform as : 5 — > 6 — AT/rj, v —>■ v — L' , 
so that two gauge invariant scalar perturbations for the radiation fluid can be constructed, 
e.g. : 

5^ = 5 + AC , v'^ = v + E'. (10) 
Instead of 5^ one can also use : 

6^ = 6--(B + v) or 6^ = 6--(B-E'). (11) 
V V 

5^, 5'^ and 5^ are the density contrasts which arc respectively defined in the comoving gauge 
where ST^ = 0, the newtonian (or longitudinal) gauge, and in the fiat slicing gauge. They 
are related by 

= S^ - -v^ and 6^ - 4*. (12) 
V 

Finally ©^j^ is the energy-momentum tensor of the scalar field at the origin of the 
topological defects. We suppose that it is a small perturbation which does not contribute 
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to the background (this is the so-called "stiff approximation" (see e.g. [8])). We decompose 
its scalar components as : 



The four source functions , ,11^ will be discussed later. They are gauge invariant 

since Q^i, is a tensor which vanishes in the unperturbed background. 

Having defined gauge invariant scalar perturbations for the metric (eq (7)), for the 
radiation fluid (eq (10-12)) and for the sources (eq (13)), we now write their evolution 
equations. Wc shall write them in Fourier space, the Fourier transform of any function 
f{x'^,r]) being defined as 



The conservation equations for the radiation fluid, when the background is governed 
by equations (2-3), can be cast under the form (see ref [5-7]) : 




(13) 




(14) 




(15a) 




(156) 



The conservation equations for read : 







(15c) 



V 





For all but the k — mode, the linearised Einstein equations read 



* - $ = kH* 



(15e) 




(15/) 




(15^) 
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In section 6 we shall solve this set of eight equations for the four unknowns ^,^,5^ 
and : the source functions p*, P^, w®, ft*, subject to the constraints (15c) and (15d) 
being known, eq (15h) will give the metric perturbation Then $ is given by (15e), 5'' 
by (15f) and i)^ by (15g). The last two equations, (15a) and (15b) are redundant (Bianchi 
identity) and can be used as a check of the calculation. 

2.3 Vector perturbations 

The line element for "vector" perturbations reads 

ds^ = alrf" [-drf + 2Bidx'dr] + {Sij + 2%^^) } dx'dx^] , (16) 

where fi* and are small functions of space and time subject to the condition : = 
diE^ = 0. (Henceforth all "barred" quantities will be divergenceless vectors : diV^ = 0.) 
Under the infinitesimal coordinate transformation 

rj^T] , x'^x' + L' (17) 

where is an arbitrary first order divergenceless vector, the four components of the two 
vector perturbations Bi and Ei transform as 

Bi ^ Bi + L[ , Ei ^ Ei + Li, (18) 

so that the two components of the vector perturbations 

^i = E[ - Bi (19) 

are gauge invariant. 

The non-zero vector components of ST^^y, the perturbation of the energy- momentum 
tensor of the radiation fiuid, can be expressed in terms of v*, the vector perturbation of 
the fiuid velocity, as (see ref [5-7]) 

K5T^i = -^{Wi + 4vi), kST^ = ^d^iEiy (20) 

In a coordinate transformation : — > — L*', so that 

vl =Vi + Bi (21) 

is a gauge invariant quantity. 

As for the non-zero vector components of the energy-momentum tensor of the scalar 
field at the origin of the topological defects we write them as 

el = -vt and ej< = 2c»(,nj) (22) 

where the four source functions vf and n| are gauge invariant and will be discussed in 
section 5. 
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The gauge invariant vector perturbations for the metric (eq (19)), for the radiation 
fluid (eq (21)) and for the sources (eq (22)) having thus being defined, we now write their 
evolution equations. 

The equation of conservation for the radiation fluid is (in Fourier space) 

= 0, (23a) 

and is 

K + - = (236) 



for the defects. The Einstein equations split into 

-A;2|. = --^^«-2avS' (23c) 

I- + = 2kuI {23d) 

The four source functions and 11^ being known and subject to the constraint (23b), 
cq (23d) and (23c) will yield the metric perturbation $j and the velocity perturbation vi. 
Eq (23a) is redundant. 

2.4 Tensor perturbations 

The tensor perturbations of the geometry are defined by 

ds'^ = alrf[-dr}^ + (% + 2Eij)dx'dx^], (24) 

with diE'^^ = El = 0. Eij is gauge invariant. The non-zero tensorial components of the 
perturbation of the energy-momentum tensor of the radiation fluid reduce to : 

^^^S = ^E,j (25) 

and the evolution of E^j is given in Fourier space by (see ref [5-7]) 

E.i + ^E[i + k^Eu^2Kul, (26) 

^ s 

where Xlj^ is the tensorial part of the energy-momentum tensor of the scalar field giving 
rise to the topological defects. 

3. The phase transition 

Topological defects are formed in the very early universe during a phase transition 
when some scalar field acquires different vacuum expectation values in different regions of 
space [2] [3]. 



Before the phase transition, the field is zero everywhere and its "false vacuum" po- 
tential energy is constant. Its energy- momentum tensor G^^, can then be modelled by a 
small cosmological constant. Hence, the only non vanishing source functions are 

p^ = Xn\ P' = -Xt,\ (27) 

where A is a positive constant related to the false vacuum energy. They satisfy the con- 
straint (15c) (cq (15d) is empty since does not depend on the space coordinates). 

Before the phase transition then, the universe is strictly homogeneous and isotropic. 
Hence the separation made here between the background governed by the radiation fluid 
and the perturbations caused by the scalar field in its false vacuum state is artificial. 
Indeed a strictly homogeneous and isotropic perturbation can always be absorbed into a 
redefinition of the background scale factor. We introduce it however for the sake of clarity. 

As the temperature of the radiation fluid drops below a given critical temperature the 
held rolls down or tunnels to a "true vacuum" state of zero energy. Bubbles of true vacuum 
are formed, grow and percolate. If the true vacuum manifold is degenerate topological 
defects appear at the boundaries of bubbles characterized by different true vacuum values 
of the held [2-3]. 

The precise energy-momentum tensor of the fleld after this phase transition depends 
on the nature of the defects formed and can be obtained solely by means of heavy numerical 
calculations (see e.g. [9-11]). However some of its statistical properties can be inferred from 
general arguments (see below Section 5). 

The duration A?] of the phase transition itself must be brief (if it is delayed for too 
long the universe may enter an inflationary phase [12] [1], and we assume that this does 
not happen). In fact we assume that it is less or of the order of one Hubble time 

A?7 ~ 1/n with n = (28) 

Vpt 

where rjpT is the epoch of the phase transition. Hence the phase transition will look 
instantaneous to perturbations evolving on time scales greater than Ar] that is to modes 
k such that 

1/A; > Ar/ 4^ krjPT < 1 (29) 

that is to modes which are larger than the horizon at the epoch of the transition. Now all 
modes of interest today were larger than the horizon in the early universe. It is therefore 
justified to describe a transition which lasts less than a Hubble time as instantaneous. 

Thus, we only need to match the spacetime geometries and the matter variables on the 
surface of transition that is the surface E of constant temperature (or constant density). 
This will be done in the next section by imposing that the induced three metric on E and 
the extrinsic curvature of E must be continuous (see e.g [13]). 

4. Matching conditions on a constant energy density surface 
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Following [14] , we first write the matching conditions in a gauge where E is a constant 
time hypersurface and then translate them in an arbitrary coordinate system. (The differ- 
ence with [14] is that here the background geometry evolves smoothly during the transition 
so that Ti and Ti' are continuous.) 

4.1 Scalar modes 



In an arbitrary coordinate system in which the line element is given by (4) the surface 
of transition S is defined by q{'q, x^) = qo{r]) + 5q = Const where q is a scalar (we shall 
take q to be the total energy density). In the coordinate transformation (5) Sq — > Sq + TqQ. 
Choosing T = Sq/qQ (and L arbitrary) therefore takes us to the gauge where the surface 
E is a constant time surface. In that gauge, that we shall label with a tilde, the unit 
normal vector to the constant time hypersurfaces is 

77-0 = — a(l + A) and = 0. (30) 

The induced metric and the extrinsic curvature are defined by 

-^^ll^= Qfiv + riij^Uy, and i^^^. = ~^^n ^,xv, (31) 
where C stands for the Lie derivative. Their non vanishing scalar components are given by 

Uj^alri\l + 2C)5,,+2d,d,E] 



1 



1 T 



V 



6)+( — A + C] 5) - d'dAB - E') 



(32) 



At first order, when one shifts back to the original, arbitrary gauge (using (6)), the 
continuity of _L^iy and K^y then imposes 



C- 



1 5q 



d'dj E' 



VPT go 

5q 



= , [a%(E + L)]^ = 



= , A + C" + ^ J 

VPT VpT % J 



(33) 



= 0, 



where [F]± is defined as [F]± = lim£^o+ [F{vpt + e) — F{vpt ~ e)] and where vpt is the 
conformal time at which the transition occurs and q is the surface of constant total density, 
so that 

g^=p' , 5q = p5 + p'/a\ (34) 

The second condition in (33) is empty since one can always choose L such that it is fulfilled. 
As for the other three, they can be rewritten in terms of scalar gauge invariant quantities 
and in Fourier space, as 



I 5^ + up' 


= 0, 


k^kj ^ 




± 





.VpT 



5* + 



Kp^ 



Vpt 



2r]PT 



+ ^kp'vIt 







= 



(35) 
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and it is an easy exercise (which makes use of the hnearised equations (15)) to show that 
the third is redundant. 

To summarize, the two independent matching conditions for the scalar perturbations 
can be simply written as (using the relation (12) between S'^ and S'^) 



+ Kp' 











- ± 





0. 



(36) 



4.2 Vector and tensor modes 

The normal vector to the constant time hypersurfaces does not have any vector nor 
tensor component and thus, 

no = —a and rii = 0. (37) 

This leads to the following expression for the vector components of the induced three 
metric and the extrinsic curvature, 



1 



As for their tensor components they are 

The matching conditions then reduce to, in Fourier space 

Nli)]± = and [4j]± = ; [^1,]± = 0. 



aQT} 



(38) 



(39) 



(40) 



4.3 The case of a phase transition 

In the particular case of the phase transition described in Section 3, the perturbations 
are strictly homogeneous and isotropic before the transition so that their Fourier trans- 
form is a zero mode proportional to the Dirac distribution 5{k). After the transition, all 
perturbations depend on space (apart from the zero mode which can be absorbed into a 
redefinition of the background). In Fourier space then the matching conditions are, for all 
modes apart the strictly A; = one : 



2 



5^ + Kp' 



= 



iir]PT 







-I r]pT 



= 



Vpt 



(41) 



ijiVFT 



0. 



There are two conditions for the scalar modes, one for the vector ones and two for the 
tensor ones, which is the required number since the equations of evolution for ^ and Eij 
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are second order and the one for $j is first order. Note that these conditions, because we 
assumed that the universe was strictly homogeneous and isotropic before the transition, 
do not depend on the physics before the transition, for example on the value of the false 
vacuum energy. 

What needs to be done now is to propagate the matching conditions (41) to a much 
later time in the radiation era, in order to set analytically the effective initial conditions 
which must be taken in the numerical integration of the evolution equations for the per- 
turbations. In order to do that we need to specify the energy-momentum tensor of the 
sources (Section 5) and analytically integrate the evolution equations from the epoch of 
the phase transition to the time when the numerical integration starts (Section 6). 

5. The energy-momentum tensor e^,. of the sources 

5.1 Coherent vs independent sources 

When the phase transition occurs the scalar field settles randomly in its true vacuum 
state in uncorrelated spatial domains. The distribution however is statistically homo- 
geneous and isotropic because the background is so and the physics of the transition is 
supposed to obey the cosmological principle. 

The ten components of the energy-momentum tensor Qp,i,{'r],x^) of the topological 
defects created are therefore ten random fields, as well as their ten Fourier transforms 
0^t/(r7, fc*) (which are complex but such that Q^^iyij], fc*) = Q^yij], ~^*))- Prom now on, we 
ignore the {k = 0) mode which can be absorbed in the background. 

The statistical properties of these ten random fields, that we shall denote collectively 
by Sa{r]., x^) or <S'a(?7, /c*), are determined by the values of their correlators. Since we have 
in view the computation of the two-point correlator of the temperature anisotropics of the 
microwave background, the quantities we need to know are the unequal time two-point 
correlators of the sources, that is 

{S,{v,x^)Sa'{v',x'^)) ^ Ca,a'{V,V',r) (42) 

where (...) means an ensemble average on a large number of different realisations of the 
transition, and where the correlator Ca,a' depends only on r = |x — x'| because of the 
homogeneity and isotropy of the distribution. The power spectrum of the correlator Ca,a' 
is defined as 

Pa,a' iv, V', k) = {27r f/^Ca,a' {v, v' , k) (43) 

where a hat denotes a Fourier transform (see (14)) and where the dependence in k = \k\, 
as well as the fact that Pa,a'{Vi v' i k) is real are again due to the homogeneity and isotropy 
of the distribution. The power spectrum is related to the correlators in Fourier space by 

{S: (ry, k')Sa' iv', k")) = 5{k' - k'')Pa,a' iv, v',k). (44) 

As shown by Turok [15], a clever way to look at the power spectra Pa,a'{Vi'n' ik) is 
to see their ensemble as a matrix M where the columns are labelled by the indices (a, 77) 
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and the rows by the indices {a',r]'). Because the power spectra are real this matrix is 
symmetric. Symmetric matrices can be diagonalised (we ignore here the fact that M is of 
infinite dimension) so that we have 

Pa,a'iV, fc) = E >^af,P%{k) (45) 

where Xaf, are the infinite number of possibly degenerate and positive eigenvalues of the 
matrix M. 

Depending on the numerical value of the eigenvalues A^^ and the associated (nor- 
malised) eigenvectors p^jj, the sources can be called coherent, independent, or mixed. 

Perfect "time coherence" of a given random field Sa-, say, means that its power spec- 
trum factorizes, that is 

Pa,aiViV'ik) =paiv,k)pair]'ik). (46) 

(This means that in (45) Xafj is of the form Xafj = X5aaS(fj — fj), and pa(r],k) = y/Xp'^^.) 
Paiv^ k) is the square root of the power spectrum of the equal time auto-correlator of Sa '■ 
(5^(77, k^)Saiv. k'^)) = 5{k^ - k'^) (paiv, k)f. 

In contrast, perfect "time independence" means that 

PaAv. rl, k) = S{rj - ri')Pa{v, k). (47) 

(This means that in (45) A^^ = AJaa, Pa'^ik) = PariS{fj - ry), and that Pa{r], k) = X {Par,f ■) 
Textures, which evolve fairly smoothly, tend to be time coherent, whereas local cosmic 
strings, which undergo complex processes of intercommutation, tend on the other hand to 
be time independent (see e.g. [16]). It is clear that the evolution of time independent 
sources is more difficult to describe since at each moment a new realization of the random 
fields is drawn. In the following we shall consider time coherent sources only. 

Time coherence (or independence), which concerns autocorrelators only, does not 
however characterize completely the statistical properties of the sources. Let us then turn 
to cross correlators. 

A given subset {Sa, a € A} of the ten sources is statistically coherent if 

Pa,a' iV, V, k) = Paiv, k)pa> [f]' , k). (48) 

(This means that in (45) the eigenvalues Agrj are zero for all a not in the subset (a, a') and 
all fj not equal to some fj.) (In keeping with the terminology used in quantum mechanics, 
such sources Sa can also be called "pure", see [15]). Note that statistical coherence of 
the random fields among themselves, that is property (48), implies time coherence, that is 
property (46). 

Knowing the unequal time two-point correlators (42) or (44), subject to the coherence 
condition (48), is not enough a priori to specify completely the random fields Sa{ri,x'^) of 
the subset. However a large class of random fields satisfying (48) is 

Sa{v,k')=Pa{v,k) eA{k') (49) 
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where eA{k'^) is a normalized complex random variable characterising the subset A : 
(e^(A;*)e^(A;'*)) = 6{k'^ — k'^). It is clear that as long as one computes two-point corre- 
lators only the two definitions (48) and (49) are equivalent. (Beware of the fact that the 
left-hand side of (49) is not some "typical" realisation of the random field but the ran- 
dom field itself, as should be clear from the explicit introduction of the random variable 
eyi(A;*) in the right-hand side. Omitting e^(/c*) on grounds of short-hand notations could 
for example induce into thinking that Saij], /c*) is a real quantity !) 

Consider now two disjoint subsets, {Sa-, a G A}, {Sb, b G B}, each being coherent, that 
is such that Sa is given by (49) and St, by a similar expression : Sb{r], /c*) = PbiV: k) esik'^), 
esik'^) being the random variable characterising the subset B, with (e^(A;*)es(A;'*)) = 
6{k' - k''). 

The two subsets A and B are said to be statistically "independent" if {e\{k'^)eB{k'^)) = 
SAB^k'-k'^). 

The difference between coherent and independent sources is most strikingly seen from 
the power spectrum of the sum of two different random fields. In the first case (a and 
a' eA) 

(K(r7, k^) + S:,{v, k%Sa{v, k'^) + 5^.(77, k'^)]) = [p«(r7, k) +Pa'{v, k)f 5{k' - k'^) (50) 
whereas in the second (a G ^, 6 G -B) 

([5:(r7, k^) + Sl{r^. k^)] [Sa{r]. k'^) + Sb{v, k'% = {\pa{v, k)f + Mv, k)f) S{k^ - k'^). (51) 

In the following we shall consider statistically independent subsets {A^ B...) of statis- 
tically (and time) coherent sources (S'a, S^, ■■■)■ Indeed, we shall see that causality imposes 
that the sources be divided into statistically independent subsets and we shall assume that 
within each subset the sources are coherent. This covers a large class of sources. Indeed, 
according to (45), the correlators of more complex, "mixed" (or partially coherent) sources 
can be decomposed into a sum of correlators of coherent sources. 

5.2 Scaling properties 

In statistical physics a dynamical system described by an order parameter 'ip{r],x'^) is 
said to follow a scaling law if the time evolution of its statistical properties depends only 
on a time-dependent length scale L{'r]), so that its equal time auto-correlation function can 
be written in the form 

(^(r7,x^)V(?7,x^)) = F(r/L(77)). (52) 

(The order parameter is supposed to be an homogeneous and isotropic random field.) In 
Fourier space the scaling law (52) translates as 

{iP*{r], k')iP{r], k")) = S{k' - k"){27rf/^L%r))F{kL{r])). (53) 

In a cosmological context a natural length scale is the Hubble radius : L{r]) = = rj. 
Dimensionless "order" parameters 'ip{r],x'^) describing the network of topological defects can 
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be constructed using background quantities, such as Q'^{r],x'^)/p, where p oc l/ry^ is the 
background energy density. Numerical simulations of texture (sec e.g. [9] [10]) and cosmic 
string network (see e.g. [11]) evolution (as well as qualitative arguments (see e.g. [7])) 
have indicated that indeed these dimensionless random fields (but not the random fields 
themselves) do obey, soon after the transition, a scaling law (in fact, if they did not they 
would be either irrelevant or catastrophic) with the Hubble radius as length scale. 
We therefore have 

{&,.{rj,x^)e,.{ri,x'^)) = ^F^Ar/v) (54) 

so that the power spectra are of the form 

(e;,(r7, k^)&^Av, k")) = S{k' - fc'O (271)3/2 lF^,(fcr7). (55) 

For statistically coherent sources described by (49) the scaling property (54-55) trans- 
lates as 

Sa{v,k') = ^Ukr])eA{k'). (56) 

What remains to be determined then is the behaviour of the ten functions fa{kr]), as well 
as the subsets A. 

5.3 Causality constraints 

The detailed structure of each realisation of Q^jy{rj,x^) is of course very complex. On 
scales less than the correlation length (of the order of the Hubble radius at the epoch of 
transition), 0^1/(^7, x*) is almost zero since the scalar field is there in a true vacuum state. 
In fact, it is non zero only around the location of the topological defects, decaying more 
or less quickly to zero away from them, depending on whether the defects are local or 
global [2]. Numerical simulations have shown that, as time passes, local cosmic strings 
interconnect and produce loops which decay by emission of gravitational radiation, in such 
a way that the energy distribution scales like the background density (see previous &). In 
the case of textures, previously uncorrelated regions become causally connected and the 
scalar field tends to take the same true vacuum value in each region, so that the number 
density of textures decreases inversely to the horizon volume ; their energy is also redshifted 
away so that, again, the energy density distribution scales like the background density. 

Now, even if each realisation of the phase transition produces a network of defects 
whose typical length scale is, at all times, greater than the horizon, their ensemble is, 
because the network appeared at a definite time, that is for causality reasons, completely 
uncorrelated on scales larger than the horizon. Therefore, as stressed e.g. by Turok [17], 
the unequal time correlators in two points P{r],x^) and P'{r]',x'^) are strictly zero if the 
past light-cones of P and P' do not intersect on the surface of the phase transition : 

{e^Ari,x')Qxp{v',x")) = if \x-x'\>ri + ri'. (57) 
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Property (57) translates in Fourier space into the fact that the power spectra are white 
noise on super horizon scales (that is for krj <^ 1). Indeed, because the correlators (57) 
have compact supports their Fourier transforms are analytic in /c*. 

Let us first consider the spatial components of the correlators. They can be written, 
in the small k limit, as 

(e*,.(r7, k')^^. k")) = 5{k' - k") [A5ij5M + B {Si^Sji + 5u5jk)] + (58a) 

where the coefficients A and B are independent of /c* (but may depend on r] and r]'). The 
right hand side term was obtained by requiring a constant tensor (with respect to A;*) which 
respects the index symmetries of the correlator (see Turok et al. [18]). More generally, the 
correlators are of the form 

(e*/??, k')@ki{'n. k")) = d{k' - r) {etU.tki + 2e^u,^kni)j) , (586) 

with 

tij = todij + tikikj, Uij = uo6ij + uikikj, (58c) 

where to, ti, uq and ui are analytic functions of k^ (e^ = ±l,e„ = ±1). This is the only 
decomposition which is compatible with the homogeneity and isotropy of the distribution. 

In Fourier space, any rank 2 symmetric tensor can be decomposed into scalar, vector 
and tensor parts, respectively : 

G,,- = eg + eg + eg (59a) 

with 

eg ^ QP.^P'^^ + L^L^^ &ki (596) 

eg^(if4 + Lfp;)efe, (59c) 

ej = (^PtP^ - lp^JP''^ Qkl, (59d) 

where the complementary projection operators and are defined as 

Pij — ^ij kikj/k , Lij — k^kj/k . (60) 

As a exercise, one can easily check that the pressure and anisotropic stress 11* defined 
in (13) are simply given in Fourier space by the expressions 

P^ = 5'^eij/3, k^Il' = ^ (P*^ - 2V^) Qij. 

It can first be shown, using the decomposition (59) with the generic term (58b-c), that 
any correlator between scalar and vector, or vector and tensor, or tensor and scalar spatial 
quantities vanishes (this in fact is a consequence of the isotropy of the distribution). 
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Let us then consider first the correlators between scalar spatial quantities. Using the 
above expressions, one finds easily 

= (A + 2S/3), (fl"*fl") = 3BA;-^ {P'*U') = 0{k^), (61) 

where l^p^*p^^ etc stand for {P^*{ri, k'^)P^{r]', /c'*)) and where it is understood that all the 
correlators are proportional to d{k^ — k'"^). 

The correlators between the other scalar terms are obtained by introducing the cor- 
relators mixing time and spatial indices, namely 

(e*oeoo) = C, (e*oei,) = DoS^j + D^k^k^, {Q*oAj) = EoS.j + E^k^kj, (62) 

(0^o©o^) = -iFki, {Qok^ij) = ikk {GoSij + Gihkj) + 2iG2k^i5j)k, (63) 

the nine functions {C, Dq, Di, Eq, Ei, F,Go, G1G2) being analytic in /c^. Performing the 
scalar, vector tensor decomposition we obtain 

{rp') = C, {rP') = Do + Diky3, (p^*n^) = -Di, {v'*v')=Eok-^+Ei, (64) 

and 

{p'*v') = F, {v'*P') = (Go + G2/3) + Gi/cVs, {v'*U') = -Gi - 2G2/k\ (65) 
Similarly, one finds for the nonvanishing vector correlators 

(nrn^') = Bk-^Pij, (Sf^J) = EoPij, = G2Pij (66) 

and for the tensor correlators 

{fCjf^ki) = B {PikPji + PiiPjk - PijPki) ■ (67) 

These results imply that, as a consequence of causality requirements, the quantities 
(56) describing the sources can be divided into three statistically independent subsets as 
defined in subsection 5.1 : a scalar subset including p'*, P'*, H'* and v'^, characterised 
by a random variable es{k^)', a disjoint, independent vector subset including v? and n|, 
characterised by e\/(/c*) ; finally a tensor subset containing Ilfj and characterised by exik"^)- 

5.4 Conservation equations 

The last general property satisfied by the sources is the conservation equations. These 
equations can constrain the correlators obtained in the previous subsection. Indeed , if 
one multiplies eq (15d) by and then takes the correlator, one finds that the leading 
term of £^0 must decay as r]~^. In the following, we shall neglect this decaying term which 
turns out to be negligible with respect to the scaling solution (see below) and thus take 
Eq (X k"^. As a consequence, 

{v^*v^) = Oik% = 0{e)P,j. (68) 
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In order to maintain statistical coherence among the scalar subset of variables, we shall 
now restrict our study to the case where the dominant k~'^ dependence of the autocorrelator 
of ri* (sec eq. (61)) vanishes (otherwise, 11* would go like k~'^ while and P* go like 
k^, which would imply, assuming coherence, that the cross-correlators of 11* with p* or 
would go like k~'^ in contradiction with (61) and (64). Thus, taking into account the next 
order terms in the analytic expansion, we shall assume 



(n**n*) = o(/cO), (n, n^.) = o{e)Pij, (% n,,) = ^(/c^) {p.^Pji + PuPju - PrjPki) ■ 

(69) 

Combining our results with the scaling property (55-56) completely determines on 
super-horizon scales the auto-correlators (48) or, equivalently for our purposes, the random 
fields (49). Introducing the scalar, vector and tensor decomposition performed in section 2, 
the random fields describing statistically coherent sources that scale, satisfy the causality 
and conservation constraints finally are : 

= ^h{kri)es{k'), P' = -^Mkr])es{k% 

(70) 

= -Vvf3{kv)es{k'), = r]'^^U{kr])es{k'), 

Vi = k^9u{kv)ev{k'), U- = krj^/%i{krj)ev{k'), = k^rj^/%j{krj)eT{k') (71) 

where the tcn^ functions fa{kr]),gai{k'r]),hij{k'r]) tend to constants on super horizon scales. 
On smaller scales {ki] ^ 1) they tend to zero (in an oscillatory fashion) at a rate which 
depends on whether the defects are local or global and which can be determined precisely 
solely by means of numerical simulations. However, since our purpose is to set initial 
conditions at a time when all relevant scales are outside the horizon, we shall not need to 
know these functions for kr] > 1. 



6. The long wavelength solution after the transition 

We shall now construct the long- wavelength solutions, that is valid as long as kr] <^ 1, 
of the equations of motion written in Fourier space for the scalar (eq (15)), the vector (eq 
(23)) and the tensor (eq (26)) perturbations, subject to the matching conditions (41) at 
V = Vpt- The source terms will be supposed to be given by Eq (70-71) where the functions 
fa{kr]), Qaiikr]), hijikr]) can be expanded in power series of kr] and tend to constants when 
/cr/ ^ : 

fa{kri) Aa a =1,2,. .A, Qai^ Bai a = 1,2, hij Qj when krj 0. (72) 

^ In fact there are six such functions only : see "Cosmic microwave background 
anisotropics seeded by coherent topological defects : a semi-analytical approach" 
by J. P. Uzan, N. Deruelle and A.Riazuelo 
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The long wavelength description of the topological defects is thus reduced to the giving of 
ten constants. 



6.1 Scalar modes 

The four scalar source terms (p*, P*, v*, n*) are subject to the two constraints (15c-d). 
In the long- wavelength limit (72), these constraints read : 

Ai = -6A2 , A3 = Ia2. (73) 

o 

The general long wavelength solution of Eq (15h) for ^ is (taking into account (73)) 

* = *o+^ + ^«?7'/'(A4 + ^2) (74) 

where and are two constants of integration and where, from now on, we omit writing 
the normalised random field es{k'^) to which each scalar perturbation is proportional to. 
Eq (15e-f-g) then give the other perturbations 



77-^ 9 5 



(75) 



Equation (12) then gives 



= -6*0 , = -2*0 + ^ + l^V^^\A2 + ^4). (76) 

(One can check that the Bianchi identities (15a-b) are identically satisfied by the solution 
(74-75).) Finally the constants and are determined by the matching conditions 
(41): 

*o = -J«^P?^2 , 4^ = J«4t(-2^4 + ^2). (77) 

One therefore sees that deep in the radiation era but long after the transition, when 
V ^ Vpt, all terms proportional to and in the expressions for ^,^,5^,v''^ and 5'^ 
can be neglected. We shall comment on 5^ in the paragraph on compensation. 

6.2 Vector and tensor modes 

The two vector source terms (vj ,11^), given by (71-72) are constrained by the conser- 
vation equation (23b). In the long wavelength limit this leads to Bu = 0. At next order 
one gets 

vt = Iv'^'k'B,,, (78) 
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where, again, the normahsed random field ey(A;*) to which each vector perturbation is 
proportional has been omitted. The general long wavelength solution of equation (23d) is 

^i = ^ + UkB2iv'^', (79) 

T}^ 9 

where $oi is a constant of integration which is determined by the matching condition (41) : 

= ~KkB2ir]p^T- (80) 

VpT " 

Finally (23c) yields 

vl = Y^oi (81) 
(so that the Bianchi identity (23a) is satisfied). 



The general solution of equation (26) for the tensorial modes, for a source given by 
(71)(72),is 



E 



kl 



Am + 



Bki 
r) ^ 63 



(82) 



where A^i and B^i are two integration constants determined by the matching conditions 
(41): 



B 



kl 



-Kk'^Cki rjpT- 



(83) 

One sees again that deep in the radiation era but long after the transition, all terms 
proportional to the integration constants can be neglected. 



6.3 Initial conditions for Cold Dark Matter 



We assumed that at the epoch of the transition the matter content of the universe 
was a radiation fluid consisting of photons, neutrinos and highly relativistic particles, some 
coupled to the photons, some not (being "WIMPS", see e.g. [19]). Thus the density con- 
trasts that we have introduced (5'', are density contrasts of the radiation fluid. Now 
as the universe expands and its temperature drops, the WIMPS become non relativistic 
and turn into cold dark matter (CDM). In this paragraph we discuss the evolution of the 
perturbations of this pressureless component of the cosmic fluid, when the universe is still 
radiation dominated. 

The background energy momentum tensor of CDM simply is : 7^;^ = PcU^Uj^ with 
uq — —a, Ui — {a being the scale factor) and where pc is the energy density. Since CDM, 
by deflnition, interacts only gravitationaly with the radiation fluid, T^i^ is conserved so 
that p'Jpc = -^/v- 

The perturbations of T^j^ can be decomposed into scalar vector and tensor components 
which can be written as 

S%l = a^PciSe + 2A) , STi = -a^pMB + Ve) , ST^^^ = -a^ p,iB, + vi) (84) 
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(the other components being zero) where 5c, f c and are the density contrast and the 
velocity components of the perturbations of the CDM fluid. In the coordinate transfor- 
mation (5) 5c ^ 5c ~ 3T/?7. As for Vc and vl they transform Hke the radiation velocities. 
Various gauge invariant perturbations can therefore be built, such as : 

Sl = Sc + 3C , vl = Vc + E' , vl' = vi + B\ (85) 

Other gauge invariant density contrasts can be introduced, e.g. 

Sl^6c-^{v + B) or Sl^6c-^{B - E'). (86) 

Note that is the density contrast in the gauge where the radiation fluid is at rest. 

Since radiation still dominates, the density contrast and velocity perturbations of the 
total fluid are almost equal to the radiation density constrast 5 and the radiation velocities 
V and t**, so that the Einstein equations (15e-h), (23c-d) and (26) remain unchanged. 
Therefore the solution given in the previous paragraphs for the metric perturbations as 
well as for the (dominant) radiation component still holds even after the WIMPS have 
become non relativistic. 

As for the evolution of the perturbations of the CDM it is governed by the conservation 
equations of its energy- momentum tensor, that is (in Fourier space) : 



vf + = 0, (88) 
where $ is given by (75). The general solution of this system is, on super-horizon scales 

« = <^o , -'c="-^-^V+^-l^{2A,-7A,W/^ (89) 

C = ^, (90) 
V 

where Sq, vq and Vq are integration constants, which contrarily to \&o and \&i are not given 
by the matching conditions at the time of the transition. Rather they should be deduced 
from a detailed analysis of the formation of the WIMPS. As for the other density contrasts 
they are given by 

^ = '^0 + ^*0 + ^^2^^73/2 , Sl = 60 + 3^0+^ + l{A2 + A^)Kr]^/\ (91) 

Once more one sees that deep in the radiation era but long after the transition, all 
terms proportional to and in the expressions for 5^, and 5| can be neglected. As 
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for the terms proportional to the unknown constants vq and Vq they are decaying modes. 
The only new constant of potential dynamical relevance is therefore Sq. 



Let us now compare the density contrasts of CDM and radiation. The isocurvature 



perturbation is the gauge invariant quantity S = Sc — jS. From (75) and (91) we have 
that it is given, at lowest order in k, by 



S = So + (92) 

We note that in 5', just like in the density contrasts 5'^ and S'"^., the terms of the type kA^?]^/^ 
are absent. We must therefore compute them at next order in k. Injecting the asymptotic 
behaviours (75) and (89) for the velocity perturbations into (15a) and (87), we obtain : 

S = So + l^o + ^KA2{k7i)W^\ (93) 

In the comoving gauge (in which the radiation fluid is at rest), the ratio of isocurvature to 
total perturbations is : 

^ T-(M^ (94) 

and does not depend on the detailed structure of the defects [as long as, of course, they 
can be represented by the random fields (70-71)]. 

7. Comments and dicussion 



In this Section we briefly compare our results to the initial conditions currently used 
in the literature. 



7.1 "Natural" initial conditions 



In [22-23], Durrer, Sakellariadou and Gangui studied the perturbations of the radi- 
ation and CDM components of the cosmic fluid seeded by the sources (70). They chose 
"natural" initial conditions, which amounts to setting = *i = in (74-76) (they con- 
sidered the scalar perturbations only). Our work justifies their choice since, as we saw, 
the homogeneous modes (in the sense of differential equations) of the evolution equations 
become negligible fairly soon after the phase transition. Their solution however differs 
from ours : indeed they only solve the conservation equations for the radiation and CDM 
fluids (that is Eq. (15 a-b) and (87)) together with the Einstein equations (15 e-f). They 
therefore ignore the conservation equations for the sources, that is Eq. (15 a-b), arguing 
that they do not hold for "incoherent" sources. Their solution hence depends on three 
constants, Ai, A3 and A4, instead of two. However it is clear that the conservation equa- 
tions for the sources (or, equivalently, the Einstein equations (15 g-h)) must be satisfied, 
whether the sources are described by ordinary functions or random fields, and that in 
[22-23] should be set equal to — ^^i. 
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7.2 Fixing the initial conditions by means of a "pseudo-tensor" 



The evolution of the perturbations in synchronous gauges where A = B = Bi = 0, 
which are used by a number of authors, in particular Pen, Spergel and Turok [9], can 
be straightforwardly obtained from the definition (10-11) (85-86) of the gauge invariant 
quantities we introduced and their explicit expressions (74-77) (89-91). For example : 

S^^yn = So + %o-^ + l'^v'^^A, and 5^^- = -1^ + %3/2^^ ^95) 
2 7 rj^ 7 

where cq is an extra constant of integration linked to the fact that synchronous gauges are 
not completely fixed. 

Now the solution (95) can of course also be obtained by a direct integration of the 
perturbation equations (15) written in the synchronous gauge. It is an easy exercice to 
show from Eq (15) and (87), and (70) (72), that the synchronous gauge density constrasts 
for the CDM and the radiation fluid satisfy, on super-horizon scales, the following equations 

4 V ri^ s/V 

the general solutions of which are 

Uyn ^ 3^syn ^ ^^^^ ^syn ^ ^^^2 + + -A2Krf'^ (97) 

4 r]'' 7 

where Dq, Di, D2 are three integration constants. After proper identification of the in- 
tegration constants one sees that (95) and (97) are indeed equivalent, up however to the 
growing mode proportional to Dirf . In order to eliminate this growing mode, a conserved 
"pseudo-tensor" r^,^, i.e. such that d^T^} = 0, was built by Turok et al. [9] [16]. On 
superhorizon scales the 00 component of this pseudo-tensor is 



ACToo = — — + ^ (98) 



and is, because r^jy is conserved, constant. Since before the phase tansition it was zero, it 
must be zero at all times : 

Too = 0. (99) 

Now, from (15a) and (10), C' = —5' /A. Using (97) one then indeed sees that the condition 
Too = is equivalent to imposing Di = 0. 

Introducing a pseudo-tensor is therefore a way to insure that the solution (97) is 
a solution of Einstein's equations. In fact condition (99) is nothing but the relativistic 
Poisson equation written in synchronous gauge. It is not a way to choose the four true 
integration constants of the general solution of Einstein's equations (15) and (87), that 
is ^^0 and ^^i, 5o and vq- These constants are fixed by Turok et al. [9] [17] by choosing 
the solution of (87) {6^^^" + 3C' = 0) to be ^ on the grounds that before the 
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transition all perturbations were zero, and by choosing "adiabatic" perturbations, that is 
imposing 5''*^"^ = 0. This amounts to choosing \I/o = 5o = 0, which contradicts (77) but, as 
wc saw, this does not really matter since the solution at late time is not sensitive to the 
matching conditions. 

7.3 "Integral constraints" and compensation 



Trashen, in [24], introduced a new way to extract from Einstein's equations integral 
equalities which relate the total matter perturbations on a Roberston- Walker background. 
They are defined in synchronous gauge and read 



{5T^ -n5T^x^)d^x= I bUSi 

JdT. 



5Tl + n5T^ 



2 ^ 



d'^x 



(100) 



B^'dSi 



where ST^!' is the perturbation of the total energy-momentum tensor of the defects and the 
cosmic fluid, where the quantities and 5^* are linear in the metric perturbations, and 
where S is a comoving three-volume, dT, being its boundary. 

These integral equalities can be interpreted as defining the energy of the perturbations 
(see [25]) and have been much invoked in the literature on cosmological perturbations 
seeded by topological defects (see e.g. [8-9] [26-27]). To understand better their role we 
shall restrict our attention to scalar perturbations and shall rewrite them, in the context 
here considered of a radiation dominated universe, as (see (4) (9) (10-11) (13)) 



Ud^x = / &dSi , / Ux'd^x = I &'dSi with W = ^ + 3k 



aE ^E VaE V 



V 



(101) 

where and C'* arc different from B^ and B^^ because integration by parts were per- 
formed. Written under the form (101) Trashen's equalities relate gauge invariant quanti- 
ties. 

As always with conservation laws in general relativity, eq (101) give some information 
about the solution of Einstein's equations without having to actually solve them. In the 
case at hand they tell the following : if the sources of the perturbations are localised 
within one horizon volume or if, as is the case with topological defects, they are randomly 
distributed on scales larger than the horizon, then the surface integrals on the right-hand 
side of (101) vanish for S larger than a horizon volume. The left-hand sides are hence also 
zero and, for scales larger than the horizon, the integrands being almost constant must also 
vanish. Now the equation W = is just Eq (15f) on super-horizon scales, that is nothing 
but another rewriting of the relativistic Poisson equation. 

Trashen's equalities were also interpreted in terms of "compensation" [8-9] [21] [27-28]. 
The idea here is that if the surface integrals in (101) vanish, then, as shown by Abbot and 
Trashen [29], the Fourier transform of U must be of order k"^. Now, as we have seen, the 
Fourier transforms of the source functions, e.g. and v^, are white noise i.e. of order 
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k^. Consequently, 5'^ must also be white noise, perfectly correlated to the sources so as to 
"compensate" them. Again this is a rephrasing of the Poisson equation. 

Finally, the fact that 6^ is much smaller than \1/ deep in the radiation era (see eq (74- 
76)) has also been interpreted by Durrer and Sakellariadou [23] in terms of "compensation" 
but it is perhaps simpler to say that this property follows from the continuity equation 
(15a). 

In conclusion, we believe that the somewhat heavy terminology used in the literature 
of "pseudo-tensor" , "integral constraints" or "compensation" is not really required if one 
straightforwardly solves Einstein's equations. 
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